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We study the statistical properties of many Brownian particles under the influence of both a spa- 
tially homogeneous driving force and a periodic potential with period £ in a two-dimensional space. 
In particular, we focus on two asymptotic cases, £i nt <C £ and £- m t 3> £, where £i n t represents the 
interaction length between two particles. We derive fluctuating hydrodynamic equations describing 
the evolution of a coarse-grained density field defined on scales much larger than £ for both the 
cases. Using the obtained equations, we calculate the equal-time correlation functions of the density 
field to the lowest order of the interaction strength. We find that the system exhibits the long-range 
correlation of the type r~ d (d = 2) for the case £- m t 3> £, while no such behavior is observed for the 
case lint <C £■ 
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I. INTRODUCTION 

To derive a description of macroscopic phenomena on 
the basis of a microscopic model is an important prob- 
lem in statistical physics. In equilibrium systems, equi- 
librium statistical mechanics provides a clear solution to 
this problem. However, when a system is out of equi- 
librium, even in a non-equilibrium steady state, no gen- 
eral framework is known except for linear response theory 
that is applied to systems close to equilibrium pj. Here, 
let us recall Boltzmann's significant research that led to 
the genesis of equilibrium statistical mechanics. He ar- 
rived at his famous formula by analyzing the simplest 
system, a dilute gas, thoroughly. Thus, in order to ap- 
proach non-equilibrium statistical mechanics, we inves- 
tigate a simple non-equilibrium system with exploring a 
relation between microscopic and macroscopic descrip- 
tions. 

Let us consider a system in which many sub- 
micrometer particles are driven by an external force in 
a solvent. The force consists of both a spatially homoge- 
neous driving force and the force generated by a spatially 
periodic potential. It has been known that such a system 
can exhibit phenomena out of local equilibrium and that 
it can be designed for experiments due to the recent de- 
velopment of the optical technology used in controlling 
and measuring particles 0, Q ■ The system also provides 
a typical example of the so-called driven diffusive system 
[iLljj- Theoretically, there exists a case that the motion 
of the particles is accurately described by a Langevin 
equation. In this system, we can investigate macroscopic 
phenomena both by a laboratory experiment and by an 
analysis on the basis of the Langevin equation that is 
regarded as a microscopic model. In particular, the sys- 
tematic calculation of macroscopic quantities from the 
microscopic model is the first step in developing a new 
framework of non-equilibrium statistical mechanics. 

Macroscopic phenomena in driven diffusive systems 
have been described phenomcnologically within a frame- 
work of fluctuating hydrodynamics 0, la, • The dynam- 
ical variable in this description is a density field, and its 



evolution equation is assumed to take the simplest form 
under the imposed physical requirements. For example, 
it is assumed that the evolution equation for the density 
field in driven diffusive systems possesses an anisotropic 
nature and no detailed balance condition. With this sim- 
ple assumption, the existence of a long-range correlation 
in driven diffusive systems was predicted even in a lin- 
ear model plIS- Furthermore, the anomalous behavior 
of the space-time correlation function of density fluctua- 
tions has been studied by analyzing a non-linear model 
for driven diffusive systems 0, ^3 . 

It is expected that the macroscopic behavior in driven 
Brownian particle systems is described by a fluctuating 
hydrodynamic equation. We then address a problem to 
quantitatively derive the form of fluctuating hydrody- 
namic equations on the basis of a microscopic model de- 
scribing the motion of the particles. If this problem is 
solved, we can calculate the correlation function of den- 
sity fluctuations by using the obtained fluctuating hydro- 
dynamic equation. Then, the calculation result is more 
quantitative than that by a fluctuating hydrodynamic 
equation assumed phenomcnologically. In general, it is 
believed that the equal-time correlation function in d- 
dimensional driven diffusive systems exhibits the power- 
law behavior of the type r~ d in the long distance regime 
and has a short-range part that deviates from the cor- 
relation determined by equilibrium statistical mechanics. 
However, recently, it has been shown that this type of 
long-r ang e correlation does not appear in some lattice 
gases [Til. Further, the short-range part of the correla- 
tion is connected with an extended thermodynamic func- 
tion in driven lattice gases [l2L Il3[. Therefore, it is im- 
portant to calculate a concrete form of the correlation 
function for driven Brownian particle systems. 

Taking these into consideration, in this paper, we de- 
rive evolution equations of a coarse-grained density field 
from a many-body Langevin equation describing the mo- 
tion of Brownian particles under a non-equilibrium con- 
dition. We first note that the many-body Langevin equa- 
tion is equivalent to a non-linear fluctuating equation for 
the density field 0] . To the latter equation, we apply a 
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system reduction method in order to describe the large- 
scale dynamics of the density field 0,0]. As a result, 
we obtain the evolution equation of the coarse-grained 
density field, and we calculate the equal-time correlation 
function. Furthermore, with regard to the calculation 
method, we extend a standard system reduction method 
so as to analyze a stochastic partial differential equation. 
Therefore, we expect that this paper contributes to the 
progress of such perturbation theory, too. 

This paper is organized as follows. In Sec. ^ we intro- 
duce a Langevin equation describing the motion of many 
Brownian particles under an external force. We then ob- 
tain a non-linear fluctuating hydrodynamic equation for 
the density field. In Sec. IIIII we develop a perturbation 
method to derive a coarse-grained fluctuating hydrody- 
namic equation and calculate the equal-time correlation 
functions of density fluctuations. Section IIVI is devoted 
to remarks on the present study. The technical details 
are presented in Appendixes. 



II. MODEL 

We study N Brownian particles suspended in a two- 
dimensional solvent of temperature T. Let Xi be the 
position of the i-th particle. Here, i is an integer satisfy- 
ing 1 < i < N. We express the a-th component of x, as 
Xi a , with a = 1,2. That is, Xi = {xn,Xi2). In addition, 
x = (xi,X2) indicates a position in the two dimensional 
space, where < Xi a < L, and periodic boundary condi- 
tions are assumed for simplicity. Each particle is driven 
by an external force fe\ = (/, 0) and is subject to a pe- 
riodic potential U{x\) with period I such that Ljl is an 
integer. For simplicity, we assume that the periodic po- 
tential is independent of x-i- Furthermore, we express the 
interaction between the i-th particle and the j-th particle 
by an interaction potential u(xi — Xj). We assume that 
the function u(r) decays to zero with a typical length, 
^int- p — N/L 2 represents the average density of the 
Brownian particles. 

The motion of the i-th Brownian particle is assumed 
to be described by a Langevin equation 



dXi, 



(It 



f 



dxi 



dUjxg) 
N 

- E 



u(Xi - Xj) +Ri a (t), (1) 



where 7 is a friction constant and Ri& is zero-mean Gaus- 
sian white noise that satisfies 



(R ia (t)R jf3 (t r )) = 2 1 T8 aP 5 l3 5{t - t') 



(2) 



It should be noted that without the periodic potential 
U, the system is equivalent to an equilibrium system in a 
moving frame with velocity f /j. Thus, the periodic po- 
tential is necessary for investigating the non-equilibrium 



nature. Note that such a periodic potential can be imple- 
mented experimentally by using an optical phase modu- 
lator 0,13. 

Now, wc define the fine-grained density field as 



jY 



Pd 



(x,t)=J25 2 (x- Xi (t)). 



This satisfies the continuity equation 

2 



dp d (x,t) 
dt 



E 



dj a (x,t) 
dx a 



(3) 



(4) 



and the expression of current j a (x, t) is derived from Eq. 
JU with Eq. J5J as follows: 



j a (x,t) = 



Pd{x,t) 

1 





( 8H 




dx a 


yS<p(x) 





I Tp d (x,t) 

1 



£,a{x, t), 



(5) 



where the real-valued functional H for a function tp(x) 
takes the form 



H((p)= J d 2 xtp(x)U(x) 

+ \ J I ' d 2 y<p(x)u{x - y)ip(y) 

+T J d 2 xtp(x)[\ogip(x) - 1], (6) 

and ^ a (x,t) in Eq. (JSJ) represents zero- mean space-time 
Gaussian noise satisfying 

^ a {x,t)^{x',t'))=28 a p8 2 {x-x')8{t-t'). (7) 

Throughout this paper, the multiplication of £ a (x,t) 
with a usual function <p(x, t) is interpreted as the 
Stratonovich rule in the space variable and the Ito rule in 
the time variable. The derivation of fluctuating hydrody- 
namic equations for the fine-grained density was reported 
in Ref. |L| . We explain a derivation method of Eq. (0J 
with Eqs. ©, ®, and in Appendix E| 

In the equilibrium case when / = 0, the form of Eq. © 
is understood by the following physical considerations. 
The first term in Eq. (J5J represents a drift caused by 
the gradient of the chemical potential that is given by 
the functional derivative of the potential i?(pd) with re- 
spect to pd- Here, the functional form of H(pd) may 
be physically interpreted by noting that the third term 
in Eq. © represents the entropy contribution of non- 
interacting particles. 

The noise term in Eq. (JSJ can be understood from 
the fact that the system with / = satisfies the detailed 
balance condition. This is verified as follows. Let At be 
a small time interval and Tr(cp — > <p') be the conditional 
probability of the density profile ip'(x) at time t + At, 
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provided the density profile is ip(x) at time t. ( Tr((y9 — > 
ip') is referred to as the transition probability from ip at 
time t to ip' at time t + At.) For the case / = 0, using 
Eqs. |@} , ©, and Q, we calculate 



Tr(p -> 9?') 
1 



A^ 6XP 



-At / d^cc 



-At / <Vx 

At 
"2T 

-0((At) 2 



7 



— d z x 



4T(p(x) 
p'{x) - 



VA~ 



ip'(x) - tp(x) 



At 



<p(x) 5H(<p) 



At 



5<p(x) 



(8) 



where A 1 is the Green function of the Laplacian opera- 
tor, and J\f is the normalization constant determined by 



V V 'Tr{y -> <(/) = 1, 



(9) 



where X>0' represents a functional measure. From Eq. 
(151) . we obtain 



Tr(y? -> y/) 
Tr((^' ip) 



exp 



gy) - g(g) 

T 



0((At) 2 



Using Eqs. © and (JTUJ, we derive 



where 



Ps(Pd) = ^exp 



H( Pd ) 



T 



(10) 



(11) 



(12) 



Eq. ifTTl) implies that Ps(pd) is a steady distribution 
functional of the density profile pd- Then, Eq. (|10fl is 
rewritten as 



Ps(ip)Tr(tp -> p') 

= Ps(</)T%' 



<^)exp(0((At) 2 )). (13) 



This is the detailed balance condition with respect to the 
distribution given in Eq. (|12fl . It should be noted that 
p ( \ in the third term in Eq. JSJ must not be replaced 
with its average value, because this replacement breaks 
the detailed balance condition. 

In the non-equilibrium case, the effect of the external 
force / is expressed only by the first term in Eq. (J5J as 
a modification of the gradient of the chemical potential. 
Although the modification from the equilibrium case is 
minimum, the modification breaks the detailed balance 
condition. Therefore, we need to analyze the evolution 
equation in order to obtain the steady probability distri- 
bution of the density field. 



III. ANALYSIS 

This section consists of four subsections. In Sec. 
IIII Al we present the basic framework of our perturbation 
method. In Sec. IIII Bl we investigate a non-interacting 
case as the simplest example. We find that there is no 
long-range correlation of the density field for the non- 
interacting systems. In Sec. IIII CI we take into account 
the effects of the particle interaction perturbatively. In 
particular, we focus on two asymptotic cases £- m t ^> £ 
and £j n t <C I. In Sec. 1111 Dl we calculate the correla- 
tion functions of the density field for the two cases and 
demonstrate that there exists the long-range correlation 
of the type l/r d only for the case ^ nt ^> I. 



A. basic framework 

We analyze the non-linear fluctuating hydrodynamic 
equation given in Eq. JI) with Eqs. (J5J, @), and Q. Our 
perturbation method is based on the expansion of the 
weak interaction, weak noise and the separation of length 
scales. In order to represent the expansion parameters 
explicitly, we replace £ a (x,t) and u(x — y) with /j,£ a (x,t) 
and Xu(x — y), respectively. We also set e = IjL. The 
parameters p, e and A are regarded as small parameters 
in our analysis. 

With this setting, we first consider the case that p, = 
and A = in the equation to be analyzed. Because 
the evolution equation is deterministic in this case, the 
density field relaxes to the steady one p^ that satisfies 

^-K^-'^sr)^ (14) 

where v s represents the average velocity of t he p article in 
the steady state, and it is expressed as EE Ei| 



Here, the function I±(xi) is defined as 



v s = -(l-e-W) 

7 Wo 



I ± (x l )= d^^^^TO-^M. (16) 
Jo 

By using the condition 

[^pTix) = p, (17) 

we derive the solution of Eq. 114|) as follows: 

pf{x)=p s {x l )p, (18) 



where 



Ps{xi) = —l-{xx). 



(19) 
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The normalization factor Z is determined by the condi- 
tion 



Jo 



dxip s (xi) = I. 



(20) 



Now, we consider the case that /j ^ and A ^ 0. Based 
on Eq. flT%)l. we set 



Pd(x,t) — p s (xi)q(x,t). 



(21) 



That is, the variable q(x, t) is equal to the average density 
p when p = and A = 0, and this variable represents the 
density modulation caused by the noise and interaction. 
Substituting Eq. J5lJ into Eqs. 10} , J5}, and © with 
the replacement explained in the first paragraph in this 
subsection, we obtain the evolution equation for q(x,t) 
as 



9q(x,t) 
dt 



-Mq(x,t) 



IK 



-A- 



1 



d 



' 3int( X : t)i 



(22) 



M=\ - 



(23) 



p s (xi) dx 

where the operator M is calculated as 

_ + T 9fog^(xi)\ 8_ 
p s (xi) 7 dxx J dxi 

T 

7 \dx\ dx\, 
and j int (x,t) is defined as 

= ^(jcOffCas.t) y d 2 y V ' p s (yi)q(y,t), (24) 

which represents a current generated by the particle in- 
teraction. 

Next, we notice the separation of length scales (e <C 1). 
We pay attention to the density fluctuations on length 
scales of order L and introduce the density field Q(X,t) 
with a large scale coordinate X = ex. While, on 
length scales of order £, the periodic potential deter- 
mines the system behavior. In order to express this in 
an explicit manner, we introduce a phase variable 9 as 
9 = mod(a;i,f) [H|l3. Obviously, U(0) = U(xi) and 

The density fluctuations with a smaller wave-number 
have a longer time scale. Then, when we focus on a 
time scale related to diffusion in the entire system, we 
assume that Q{X, t) obeys an autonomous equation and 
that q(x, t) — Q(X, t) depends on time t only through the 
density field Q(X,t). These assumptions are expressed 
as 



q(x,t) = Q(X,t)+p{6,Q), 



(25) 
(26) 



where fl(Q) represents a map providing a function of 
(X,t) for the density field Q, and p(9, Q) represents a 
similar map for each 8 (see Ref. [THl 

Hereinafter, we treat 9 and X as independent variables. 
Then, when the spatial derivative acts on a function of 
(9, X), it should be read as 











dx 89 1 



d 



(27) 



Furthermore, the two quantities £(x, i) and u(x) in Eq. 
(J2H with Eq. can be expressed as £(6,X,t) and 

u(9, X), respectively, by using a method presented in 
Appendix Using this expression, Eq. Q leads to 

(U0,x,t)^(e',x',t')) 

= 25 a0 £6(6 - 9')e 2 S 2 {X - X')8(t - t'), (28) 

as explained in Appendix IB1 Further, using u{9, X), 
j int (x,t) in Eq. $Z%fy is expressed as 

=p s (9) [Q(X,t)+p(9,Q)} 
d9 ei + dx)J £ 



d 2 V 



e',x-Y) 

(29) 



Ps (e')[Q(Y,t) + P (e',Q)], 

B9|l . Finally, M can be ex- 



where we have used Eq 
panded as 



M = M (0) + eAf (1) + e 2 M (2) , 
where M^> with i = 0, 1 and 2 are calculated as 

T d\og Ps {9) \ d_ T d 2 
7 d9 J 89 + 7 d9 2 ' 



MW = 
M< 2 > = 



p s (0) 

Vs 

'pM 

T f d 2 
7 U^i 2 



Td\og Ps {9) , Td 



+ 



7 
d 2 

oxi 



d9 



d 



7 89 J dX 1 



(30) 

(31) 
(32) 

(33) 



Using these, the substitution of Eqs. an d lt2lfl) 

into Eq. yields 



n(Q) 



d \ 
d9 ei + e dx 




y/Ps(0)(Q + p(9,Q))l(9, X, t) 



-A- 



9 



ax y 



^Uint(^,t), (34) 



Ps(9) \89 

where (5p/SQ) is the operator on a function ip(x) defined 



P (9,Q + cp)-p(9,Q) 



6p_ 

SQ 



^0r) + O(M 2 ), (35) 
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in the limit \(p\ — > with an appropriate norm | - | of the 
function space. Recall that p is a map in the function 
space for each 9. Thus, the operator (Sp/5Q) mathe- 
matically corresponds to a Frechet derivative when the 
function space is properly defined. This should not be 
confused with the functional derivative used in Eq. 

Now, we derive £l(Q) and p(9, Q) by solving Eq. 
with the perturbation method. Because we assumed that 
e, p, and A are small, we expand p(9, Q) and £l(Q) in 
these parameters. Specifically, setting p — e, we assume 
the form 



p(9,Q) = e Pl (9,Q) + e 2 p 2 (9,Q) 
fi(Q) = eil^Q) + e 2 n 2 {Q) + ■ 



(36) 
(37) 



Furthermore, for both pi(9,Q) and fli(Q), we consider 
an expansion with regard to A. In summary, we analyze 
Eq. (EU) with Eqs. ®, JU, and p = e and 

derive pi(9,Q) and £li(Q) iteratively. These calculations 
are shown in Sees. IIIIBI and lHI CI 

In order to make our perturbative calculation concrete, 
we introduce a space T consisting of all complex- valued, 
square-integrable, periodic functions of 9 on the interval 
[0,1]. We endow this space with the inner product 



(a, j3) 



A9 



a* (6)13(6) 



(38) 



for a,/3 e f, where * denotes the complex conjugation. 
Since the operator Ml ' is a linear map from T to T, we 
can define all the eigenvalues Xj and the corresponding 
eigenfunctions $j(&) of the operator in T by the 

equation 

MW<t> j (9) = \ j $ j (6), (39) 

where the index j = 0, ±1, ±2, ■ • • is determined by the 
relations Xj = X*_j and A = > Re(A±i) > Re(A± 2 ) > 
■ ■ • . When a complex eigenvalue is degenerate, the cor- 
responding labeling is modified appropriately. Because 
Af(°) is not a Hermitian operator, it is convenient to in- 
troduce the adjoint operator of M*- ^ through the rela- 
tion 

(M^ta,/?) = (a,Af( 0) /3) . (40) 

In this space T ', is explicitly represented as 

d_ ( v s _ T d\og Ps (9) Td_ 
89 \p s (9) 7 d9 7 d9 / 



M(°)t 



(41) 



Note that the set of eigenvalues of M^ ^ is identical to 
that of M<°). Then, we denote the eigenfunctions of 
jVf(°)t as ^j(9) and choose their labeling such that 



M {Q) ^ j (9) = A**j(6»). 



(42) 

We can choose these eigenfunctions such that the follow- 
ing hold: 



(*i,*i)=* 



(9)^(9') 



--6(9 



(43) 
(44) 



In particular, we choose the zero eigenfunctions as 

*o(0) =*.(*), (45) 
*„(0) = 1. (46) 

B. non-interacting system 

We investigate the non-interacting system by setting 
A = in Eq. J2H) with Eq. Equation J2H) can then 

be rewritten as 

SQ 

M[Q + P (9,Q)} 



n(Q) 



~p s (6) ]l 7 {dd Gl+e dX 

(VpM [Q + p(9,Q)W,x,tj) , 



(47) 



for which we use the expansions given by Eqs. I|3UI) . 13611 . 
and (|3T|) . Selecting all the terms proportional to e in Eq. 
(P7Jl . we obtain 

fii(Q) =mW Pi (6,Q)+mWQ 



(48) 

Here, f2i(Q) and pi(9,Q) are unknown functions. In or- 
der to derive them, we set 

s%(0,Q) 



+ 



^^|(v^M^)a(^,t)), (49) 

and rewrite Eq. 148fl as 

MW Pl (6,Q) = Sl (9,Q). (50) 

We can regard Eq. H50[) as a linear equation with respect 
to pi because s\(9,Q) does not contain p\. 

bincc M(°) has a zero eigenvalue, M' ) is not invertible. 
In this case there is no unique solution of pi to Eq. I|50|) 
but either no solution or an infinite number of solutions. 
Then, in order to perform the perturbative calculation 
consistently, we impose the solvability condition 

(*o,«i)=0; (51) 

under this condition, there exist solutions with an arbi- 
trary constant. This condition determines £li(Q) as 

fii(Q) 




y/ Pa Q(X,t)i 1 (;X,t) 



(52) 
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where d represents the partial derivative with respect to 
8. We note that the second term in the right-hand side 
of Eq. (|52|l is equal to zero (see Appendix [5J) • Then, 
substituting Eq. (|32|) into Eq. (|52|l . we obtain 



fii(Q) 



dQ 



(53) 



Under the solvability condition given by Eq. Ij51|l . we 
can derive the following solutions of the linear equation 
expressed by Eq. (|5U|l : 

Pi(0,Q) 

= 2.^ r (*n,-- + -a(l0g Ps ) 



dX 1 



^ifi m(* n i a[ ^ l( ., x , t)] 

•X*o(0). 
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(54) 



Here, x is an arbitrary constant. We set x = hereafter. 

Next, we will determine £l 2 (Q) and P2(#, Q)- Using the 
terms proportional to e 2 in Eq. (|47|) . we obtain 



dpi 

M ( °V 2 + M (1) pi +M (2) ( 



l 



<T d 



Ps (9) V 7 



(55) 



In the same manner as that in the first order calcula- 
tion, we impose the solvability condition for the linear 
equation of p 2 . This yields 



n 2 (Q) 

= (*o,A^ (1) Pi) + (*o,M (2) Q) 

1 



IT d 
~dX 



y/Ps 



(56) 



Here, we calculate D, ((X,t), and rj a (X,t) as 



D 



b, 



Vs , an \ 
h —d(logp s ) 

Ps 7 



T 
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C(x,t)=J^d^,^p-M;X,t) 

Va{X, t) = J-^ (*0, -ifaO. X *)) . 



where b(6) is defined as 
b(9) 

Hi 

P: 



(58) 
(59) 

(60) 



E(*o, 

m=i0 



T 2T 
—d(\ogp B ) H d 

7 7 . 



*m(0) 



(61) 

We find that b(6) is a real function (see Eq. (jCllf! *). 

Now, we write the coarse-grained hydrodynamic equa- 
tion by defining 



Q(x,t) = Q(X,t), 



(62) 



1 



H 1 (x,t) = -[C(X,f)+r ?1 (X,t)], (63) 
e 

HaOM^i^X,*). (64) 
e 

Using these and from Eqs. I|25l) . 15311 . and (|57|l . we obtain 



dQ 

dt ~ ^ 



dJ a (x, t) 
dx a 



(65) 



with 



dQ 



Jx(x,t) = v s Q(x,t) - D— + \l Q(x,t)E 1 (x,t), 



T d(q) 

J 2 {x,t) = h \jQ(x,t)Z 2 (x,t), 

7 OX2 



(66) 



where ^ a with a = 1, 2 satisfies 

(H^sc, t^sc', 0) = 2B a/J 5 a (x - ~ 0- ( 67 ) 

The noise intensities B a p are calculated as B 12 = B 2 i = 
and 



T f dd . . 

-B22 =— • 



d f b(6) 
d0 VpTM 



1 2 



1 



(68) 
(69) 



Using Eqs. (|3~2"|) and Ip^ . we obtain 

a 2 Q t a 2 Q 



n 2 (Q) =£> 



ax 2 7 sa 2 2 

jL-^/Q[((X,t)+ m (X,t)] 
d 



ox 2 



Qm(x,t). 



(57) 



See Appendix El for the calculation. 

Here, we present two remarks on the coarse-grained 
hydrodynamic equation given by Eq. (BSJ with Eqs. JSJ 
and (|67|l . The first remark is on the expression of D given 
by Eq. I|58|) . Although the expression is complicated, we 
can rewrite Eq. (|58|l as 



D 



7 Wo 



66 



1.(9) 



^[J_(0)] 2 /+(<?) (70) 
1 
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using Eq. (|16f) . The derivation is presented in Appendix 
\C\ This expression of the diffusion constant coincides 
with that of the diffusion constant of a Brownian particle 
in the tilted periodic potential \PA ll8L IT9I ] . Physically, 
this coincidence is obvious, because we consider the non- 
interacting particles in this subsection. 
The second remark is on a special relation 



Bn = D. 



The proof of this relation is presented in Appendix [U] 
This relation corresponds to the fluctuation-dissipation 
relation of the second kind in this fluctuating hydrody- 
namic equation. Using this property, in the same manner 
as that in Sec. [HJ we can prove the detailed balance con- 
dition of the system in the moving frame with velocity 
v s . From this condition, we find the steady probability 
distribution functional of the coarse-grained density field 
Ps(Q) as 



Ps(Q) = y exp (- J d 2 xQ(x) pogQ(x) - 1 



(72) 



This implies that in the non-interacting system, the den- 
sity field does not exhibit a long-range spatial correlation 
even if the system is out of equilibrium. 



C. weakly interacting system 

In this subsection, we extend the analysis in the previ- 
ous subsection to a system consisting of interacting par- 
ticles. Concretely, we expand Pi(9, Q) and Qi(Q) in Eqs. 
(EU) and 123) as 

Pi(9,Q) =Pio(0,Q) + \pn(9,Q) + X 2 Pt2 (9,Q) + ■ ■ • , 



fii(Q) =fiio(Q) + Afiii(Q) + A 2 ^ 2 (Q) + 



(73) 
(74) 



Further, in order to make the calculation results explicit, 
we consider two asymptotic cases: (i) £i n t ^> I and (ii) 
Ant < i- By substituting Eqs. (|73J) and J73I) into Eq. 
(|34fl with Eq. (|24|) for the two cases, we determine 
Pik{9,Q) and flik iteratively. 



1. case (i) 

We study the case £i nt ^> I. Specifically, we assume 
that £i„t ~ 0(e~ 1 £), and then we set 



,X) = e 2 u L (X). 



(75) 



The factor e 2 is introduced in order to develop a system- 
atic perturbation method. Substituting Eq. I|75l) into Eq. 
(|29[1 , we express the current generated by this interaction 



potential as 



1 



7 



p s (9) [Q(X,t)+p(9,Q)] 



2v dujXX-Y] 
dX 



d 2 Y 



Q(Y,t). 



(76) 



(71) Then, Eq. (|33} with Eq. becomes 



Q(Q) + ^-Q(Q)=M[Q + p(0,Q)] 



Pb(0) y 7 \d0 dX_ 
{VPs(0) [Q + p(0,QM0,X,t)} 

A _J_ fd_ \ 

' 7 £ P S W {d9 ei+€ dXj 
Ps (9) [Q + p(9,Q)] J d 2 Y dMX_Y) 



Q(Y,t) 



(77) 



We now analyze Eq. 177(1 with the expansions given by 
Eqs. (2JJ, ESI, E3, G3, and Gil- 
Selecting all the terms proportional to e, we obtain 



fii(Q) 



1 IT d 



p s (9) V 7 09 
A 1 d Ps {9) 
'ip s {9) 09 



VM9jQ^i(9,X,t) 
Q(X,t) 



^Y dUL f x ; Y) Q(Y,t). (78) 

The terms independent of A reproduce Eq. Thus, 
^io(Q) and pio are equal to fli(Q) and pi given by Eqs. 
(|53[) and l|54|) . respectively. Next, extracting all the terms 
proportional to A in Eq. (|78l) . we obtain 



m^pu= nn(Q) 

1 1 0p B (0) 
7ft(#) 50 

2 8u L (X-Y 



d z Y- 



Q(X,t) 

Q(Y,t). (79) 



dXx 

The solvability condition for the linear equation of p\\ 
yields 



O u (Q)=0. 
Then, we can solve pn as 



(80) 



Pi 



1 {0,Q) =-Y,^-(^n,ldp )QiX.n 



7 ^ A « 



d3y s« L (jr-y; 



(81) 
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where the term proportional to &o(0) is set to zero. 

Next, using all the terms proportional to e 2 in Eq. (|77|) 
with the expansions, we calculate £2 2 o(Q) and ^21 (Q) by 
repeating the same analysis. Obviously, f22o(Q) is equal 
to VL 2 {Q) in Eq. I|57[) . Extracting the term proportional 
to Ae 2 in Eq. O, we obtain 



=M (0) P2i +M (1 Vn 




V 7 50 

1 1 a 



1 d 



\p a (d)pio(0,Q)] 



dX x 



■Q(Y,t) 



Q(X,t)Jd*Y du ^f x Y) Q(Y,t) 



(82) 



Then, the solvability condition for the linear equation of 
P21 yields 

nai(Q) 

= (*o,M (1) Pn) 



7 ax 



(83) 



Substituting Eq. I|81|) into the first term in the right-hand 
side of Eq. leads to 



Pn 



=- (b, —dp, 
7 V Ps 



d 



Q { X,t) I - Y ' 



Q(Y,t) 



(84) 



where b(8) is given by Eq. Here, we note the iden- 

tity 




(85) 



where we have used Eq. (jClfl to obtain the third line, 
and the fourth line can be confirmed directly from Eq. 



(HSJ) (see also Ref. ,19j). Using Eqs. (JHU and (JHSJ, we 
rewrite Eq. (|83|l as 



n 2 i(Q) 

dv= d 



df dX x 

1 a 



7 dX 2 



Q(X lt )J^Y d ^f x - Y) Q(Y,t) 
Q(X,t)Jd*Y ^f x - Y) Q{Y,t) 



Finally, we derive the coarse-grained hydrodynamic 
equation from Eqs. (gSJ, JHOJ, and Us- 

ing the variable defined in Eqs. I|62|). i|t)3[l . and l|64|l. we 
obtain the continuity equation for Q(x,t) expressed by 
Eq. (|65|l with the current as follows: 



Ji(x,t)=v s Q(x,t)-D^-+jQ(x, t)S x (as, t) 
+ A^(l-*)Q(x,t) | d 2 y^-^g(y,t), 

J 2 (as,<) - + \/Q(x,t)E 2 (x,t) 

7 0x2 v 



+ A^,t)/d 2 / u(a; - y) - 



1 

-1 

7 



0(»,t), 



(87) 



where 6 is the dimensionless parameter defined as 
, , T dv s 



The Einstein relation leads to 6 = when / = 0, while 
S 7^ when / 7^ as far as we checked numerically (see 
Ref. H3). 

When S = 0, by using a similar argument as that in 
Sec. m we can prove the detailed balance condition of the 
system in the moving frame with velocity v s . However, 
when 5^0, the coarse-grained hydrodynamic equation 
does not possess the detailed balance property for any 
moving frame because we cannot construct a potential 
function such as H for the argument in Sec. [H] It should 
be noted that the noise intensities are not modified by 
the lowest order contribution of the interaction effects. 
Therefore, the fluctuation-dissipation relation of the sec- 
ond kind is maintained in this hydrodynamic equation. 



2. case (ii) 

Next, we study the case t mt i- We assume the form 

u(6, X) = u a e6(9)e 2 5 2 (X), (89) 

where the intensity uq is determined from the interaction 
potential with £- mt -C I. Here, we set uq = eus in order to 
develop a systematic perturbation. Note that the quan- 
tity uq should appear when we calculate experimentally 
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measurable quantities. Substituting Eq. I|89|l into Eq. 
(|29[) . we obtain 

=e^p s (6) [Q + p(6,Q)} 
1 



Then, Eq. |j3JJ with Eq. (25J) becomes 
n(Q) + ^-V(Q) = M[Q + p(9,Q)] 



(90) 



IT 1 



-Mi 



(9 



7Ps (e) \de ei + e dx) 
VpMW+WMW, x, t) 

Xu s 1 (d_ d \ 

~M9) \d~e ei + e dx) 



\p s {0) [Q + p(8,Q]\ 



de 



(91) 



We analyze Eq. I|91l) using the expansions given in Eqs. 
CP, (H3), (23, G3>, and d- 

The calculation procedures hereafter are the same as 
that for case (i). Thus, without repeating the calcula- 
tions, we summarize the results as follows: 



n n (Q) = o, 



7 t— 1 A 



de' m 



o £ p a (0') dff 
n 21 (Q)=-^-[DQ 2 (X,t)} 
where v is defined as 



(92) 

(93) 
(94) 



Tft (tf)/ + (fl)_. 



(95) 

Note that in the equilibrium case (/ = 0), v = because 
p s (6)I + (6) is equal to unity when / = 0. 

Finally, we derive the coarse-grained hydrodynamic 
equation from Eqs. (53), J^ZJ, (EH, and Using 
the variable defined by Eqs. I|u2|). . and 164|) . we ob- 
tain the continuity equation for Q{x,t) expressed by Eq. 
(|65fl with the current as follows: 

dQ 



Ji{x, t) = v s Q(x, i) - D-^- + \] Q(x, <)Hi(x, t) 



+ \vQ 2 {x,t), 

T dQ(x,t) 
7 dx 2 



+ jQ(x,t)E 2 (x,t). (96) 



Here, we have defined v as 



u _ 

V = — i/, 



(97) 



where the quantity v is independent of e and its value can 
be determined from the Langevin model we study. Note 
that the coarse-grained hydrodynamic equation does not 
possess the detailed balance property for any moving 
frame except for the case where v = 0. 



D. correlation function 

In this subsection, we calculate the equal-time correla- 
tion function Co(r) of the coarse-grained density field for 
the two above-mentioned cases, case (i) £[ nt ^> £ and case 
(ii) £i n t <C £■ We first define the space-time correlation 
function as 



C(r, t) = (ip(x, t)tp(x + r, t + r)> , 



where ip(x,t) represents the deviation of the density field 
from the average value p as follows: 



i>(x,t) = Q(x,t) - p. 



(99) 



In Eq. (|98f) . we assume that the statistical properties of 
ip(x,t) are translational invariant in the space and time 
directions. 

Hereinafter, we denote the Fourier-transformation of a 
function f(x,t) as 



f(k,u) = / d 2 xdtf(x,t)e 



-ih-x — iut 



(100) 



We also use the same notation for the Fourier- 
transformation of a function f(x) as follows: 



f(k) = / d 2 xf(x)e 



—ih-i 



(101) 



Then, we can derive the relation 

C(k,u)(2ir) a 6 2 (k + k')S(w + lu') = (4>(k,u)ip(k' ,u'Y 

(102) 

Further, the equal-time correlation function C'o(r) is 
given by 

d 2 k 



C o(r) = / 7^C (k)e ik - r 

-C{k,tu)e lk - r . (103) 



(2tt)2' 
d 2 k doj 



(2tt) 2 2tt 



In the argument below, we first calculate C{k 1 uj) from 
the coarse-grained hydrodynamic equations for cases (i) 
and (ii), and we derive Coir) using Eq. (|103fl . 
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1. case (i) 

Let us consider the case i- mt ^> I. Substituting Eq. 
(l9~9")l into the continuity equation expressed by Eq. lj6"B*|l 
with Eq. (|87|l . we can obtain the equation for ip. We fur- 
ther linearize the obtained equation. Then, the resultant 
equation becomes 



9V 9 J. 



01 



dx a ' 



(104) 



with the current J a expressed as 



+ VpHi(a;,t), 

7*2 7 J 9x 2 
+ ^/pE 2 (x,t). (105) 

The Fourier-transform of the evolution equation is 
written as 



ip(k, w) = G(k, uS) i-iy/pk ■ E(fc, oj)j . (106) 
Here, G(k,u>) is the Green function calculated as 

n! } ; =t(w + w s fci) + 9s(k)Dkj + 9o(k)-kl (107) 
where we have defined <7a(fe) as 



1 + Apu(fc) 



ffJ (fc) = 1 + _ J). 



(108) 



Using the relation 



/s a (fe,w)H^(fe',wO) = (27r) 3 2B Q/3 ^(fe + fc')(5(w + ^'), 

(109) 



we obtain 



C(fc,w) = 2p|G(fc,w)r -D/s? + -hi 

7 



(110) 



Integrating Eq. (|110fl with Eq. (|1(JT|) over the frequency, 
we calculate 



Cb(fc) =P 



gs{k)g {k) 
h A^ (fc) 1 



2T gs(k)g (k) 
g s [k)Dkl-g {k){T/i)k\ 
g 5 {k)Dkl + g {k){Th)kl 



(in: 



The asymptotic behavior in the range |fe| <C ^ in J in Eq. 
(Illl|l is evaluated as 



Co(fe) ^ 



gi/2(0) 

w(o)so(o) 

A^ (0) 



1 



2T <? 5 (0). 9o (0) 
^(0)Jfc 2 -. 9 o(0)(r/ 7 )fc 2 
9s (0)Dkf + g (0)(T/ 7 )kl 



(112) 



From this expression, the asymptotic form of Co(r) in 
the range \r\ ^> ^ nt is derived as 



Co(r) 



Ap 2 u(0)J 



27rT A /[« ?5 (0) 5 o(0)] 3 OT/ 7 

{^(o)z?]-V? + M0)(t/ 7 )]-S 2 } 2 ' 

(113) 

where r = (ri,r 2 ). This represents a long-range correla- 
tion of the type 1/r 2 . 



2. case (ii) 

Next, let us consider the case t ln t <C ^. In this case, 
we obtain the continuity equation for ip with the current 
J a expressed as 

Ji(x,t) =(v s + 2p\v)ij) - D—i- + yjp + ip(x, t)3i(!B, t) 

OX i 

+ \vi/j 2 (x,t) 



T dip 



Ja{x,t) =- —^- + y/p + i>(x,t)Ba(sc,t). (114) 
7 ox 2 

Further, for simplicity, we assume that the ip dependence 
of the noise term can be neglected. Then, the evolution 
equation coincides with the special case of that investi- 
gated in Ref. ^j|> where the fluctuation-dissipation re- 
lation of the second kind is satisfied. According to the 
result of Ref. ,10], the equal-time correlation function 
is not modified by the interaction effects up to the sec- 
ond order of A in this case (see Eq. (B2) in Ref. [lojp. 
Thus, we conclude that there is no long-range correlation 
within this approximation. 



IV. CONCLUDING REMARK 

The main achievement of this study is the derivation 
of the coarse-grained fluctuating hydrodynamic equation 
for the driven many-body Langevin system. In the two 
asymptotic cases for the interaction range between par- 
ticles, which are given by Eqs. (|75|l and ^89p. we derived 
the two expressions of particle current, Eqs. I|87|l and 
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(|96|l . respectively, with the continuity equation expressed 
by Eq. |JBS}. Using the obtained evolution equations, 
we calculated the equal-time correlation function of the 
coarse-grained density field for each case. We found that 
this correlation function exhibits the long-range correla- 
tion of the type r~ d in the case given by Eq. Ij75jl. while 
no such behavior was observed in the case given by Eq. 

We derived the coarse-grained fluctuating hydrody- 
namic equation by applying a singular perturbation 
method to a stochastic partial differential equation. The 
method is standard except for the treatment of space- 
time noise (see Appendix |B|) . We expect that our method 
can be used to investigate other related problems such as 
phase diffusion behavior in periodic pattern formations 
under the influence of noise. We also note that the de- 
rived coarse-grained fluctuating hydrodynamic equations 
in this study contain non- linear functions of / such as v s , 
D, B a , S, and v. Therefore, we can discuss the system 
behavior in a non- linear range with respect to /. 

The long-range correlation we obtained for the case 
given by Eq. (I75[l has the essentially same mechanism as 
that of Ref. [2(| which studied the system consisting of 
two Brownian particles under an external force. Thus, 
the result in this paper is regarded as an extension of 
that in Ref. [2(|, although the Fokker-Plank equation 
was analyzed in this reference. On the other hand, the 
long distance behavior for the case given by Eq. I|89|) 
might be strange, because it has been speculated that 
an anisotropic system without the detailed balance con- 
dition generically exhibits the long-range correlation 
Here, it is noteworthy that the long-range correlation of 
the type r~ d does not appear in driven lattice gases with 
the evolution rule called an exponential method, while 
it appears in the cases of a heat bath method and a 
Metropolis method . It might be interesting to find a 
connection of our result with that reported in Ref. [Tl| . 

Our calculation result for the correlation function given 
by Eq. I|lll|) provides the functional form of its short- 
range part. In contrast to the long-range part, the short- 
range part depends on the details of the system such as 
the selection of the interaction potential. Such a non- 
universal part has never been investigated intensively. 
Here, let us recall that the statistical properties of den- 
sity fluctuations are described by the free energy function 
of the system if the system is in equilibrium. There- 
fore, it might be expected that the short-range part of 
the correlation function is related to a thermodynamic 
function extended to non-equilibrium steady states. Al- 
though thermodynamics in non-equilibrium steady states 
has not yet been established, there exists one promising 
approach to construct a consistent framework whose va- 
lidity can be checked experimentally |13j. According to 
this framework, the statistical properties of density fluc- 
tuations can be described by an extended free energy de- 
termined operationally when the effect of the long-range 
correlation is removed. Indeed, by numerical experiments 
on a driven lattice gas, it was demonstrated that the in- 



tensity of density fluctuations of a particular type is de- 
termined by an extended free energy |l2j . and this was 
proved in Ref. 13]. We expect that a similar analysis 
can be performed for the Langevin system under investi- 
gation in this study. The calculation of the short-range 
part of the correlation function is indispensable in this 
analysis. 

The most ambitious goal is to provide a unified descrip- 
tion of density fluctuations in an elegant manner. Even 
if the short-range part of the statistical property of den- 
sity fluctuations is determined by an extended thermody- 
namic function, the long-range correlation that destroys 
the extensive nature of the system is obviously out of 
thermodynamic consideration. Here, it should be noted 
that a variational principle referred to as the additivity 
principle is effective to describe the long-range behavior 
of density fluctuations in non-equilibrium lattice gases 
[2ll E3 | . Although we do not know a class of models to 
which this principle can be applied, it is interesting to 
determine whether this principle can be applied to the 
Langevin system under investigation. 

The authors acknowledge H. Tasaki and K. Hayashi 
for useful discussions and comments. This work was sup- 
ported by a grant from the Ministry of Education, Sci- 
ence, Sports and Culture of Japan (No. 16540337). 



APPENDIX A: FLUCTUATING 
HYDRODYNAMICS FOR FINE-GRAINED 
DENSITY FIELD 



We derive the evolution equation for a fine-grained 
density field from the Langevin equation given by Eq. 

Mathematically, the essence of the derivation is in 
employing the Ito formula for arbitrary functions of the 
fluctuating variable 23], and the evolution equation for 
the fine-grained density field is obtained by a standard 
treatment of Dirac's delta function. This approach was 
performed in Ref. ■ Although Ref. 0] provides suffi- 
cient information on the derivation of the evolution equa- 
tion for a fine-grained density field, we present a different 
method for the derivation in this Appendix. We find that 
this method is less mathematical, but more pedagogical 
than the standard one. 

Let At be a sufficiently small time interval, and set 
t n = nAt. Then, from Eq. the movement of each 
particle during the time interval, Axi a (t n ) = Xi a {t n+ i) — 
(t n ), can be expressed as 



jAx ia (t n ) 



E 



Ou(Xi — Xj) 



dxn J 4rf dxi, 

W ia (t n ) + O ((At) 3 / 2 ) 



At 

[x e =x e (t n )] 

(Al) 
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where 



APPENDIX B: FUNCTION OF <9,X) 



W ia (t n ) = 



dtRiait), 



(A2) 



and it should be noted that the equality 

W ia {t n )W ]P {t m ) = 2 1 TS iJ S a0 S nm At (A3) 

holds almost surely |23j. 

Next, for the fine-grained density field p d (x,t) defined 
by Eq. ©, we obtain 

7 [p d (x,t n+ i) - pd(x,t n )] 
d 

= y]AXia(t n )- p d (x,t n ) 

+ \ ^2^2 AX ' Q (* n ) Ax lP( t n) JT~ Pd( x , *n ) 

+ O ((At) 3 / 2 ) . (A4) 

Substituting Eq. i|Al|l into the above expression and us- 
ing Eq. 1A3(1 . we derive 



7 [pd{x, t n+ x) - pd(x, t n )] 



= -At 



d 

dxi 



dU{xx, , 

/ KZ ) Pd{X,t) 



dx\ 



At— • / d 2 yp d (x,t) 



du(x — y) 
dx 



Pd(y,t) 



& 



d 2 



9a;j eta 



5_ 
3a; 



W 



( X ,t n )+o({Atfi 2 y 



where we have defined 



N 



W a (x,t n ) = J2w ia (t n )5 2 {x-Xi). 



(A5) 



(A6) 



Note that W a (x,t n ) satisfies 

(W a (x,t n )Wp(x',t m )) 

= 2jTp d {x, t n )S aP S 2 {x - x')S mn At. 

Finally, taking the limit At — > 0, we obtain 



(A7) 



dpd{x,t) 1 d 
7 dx\ 



Of 



dU( Xl ) 
dxi 



+ f Pd(x,t) 



~-.J d yp d{x ,t) 



du(x — y) 
dx 



Pd(y,t) 



T f(P_ cf_ 

7 \9x 2 dx 2 



2 ' aj ; Mx,t) 



S_ T Pd (x,t) 
dx y 7 



(A8) 



where £ Q (a:,i) satisfies Eq. JJJ. ft is easily confirmed 
that the final expression is equivalent to Eq. Q with 
Eq. ©. 



In Sec. IIII Al we introduced the functions £(9, X,t) 
and u(8,X) corresponding to £(x,t) and u{x), respec- 
tively. In this Appendix, we present a method to con- 
struct the function of (9, X). For simplicity, we consider 
functions defined in a one-dimensional interval, but the 
argument presented below can be extended to functions 
in two- or higher-dimensional regions. 

Concretely, we construct a function cf>(6, X) corre- 
sponding to a function 4>(x), where <fi(x) is defined in 
the interval [0, L] and <fi(9, X) is defined in the region 
[0,£) x [0,£]. We assume that there exists an integer Ni 
satisfying L — N\£. We explain a numerical computation 
method to obtain the function (f>(9, X) from <j){x) without 
a rigorous mathematical argument. 

We first divide the interval [0, L] into small segments 
[ia, (i + l)a), where < i < N a Ni - 1 and a = i/N a . For 
the function <fr(x), we set 



<j>i = 4>{ia). 



(Bl) 



Then, <pi is regarded as a real-valued function from inte- 
gers < i < N Ni — 1. For each i, there exists a unique 
pair of the integers io and i\ that satisfy 



i = ixN + io, 



(B2) 



where < i < N - 1 and < i\ < Nx - 1. Here- 
inafter, io and i\ are regarded as functions of i. Using 
this notation, we define the function <pio,ii as 

4>i M = <Pi- (B3) 

We expect that the function 4>(x) appearing in physics 
is well approximated by using & with a sufficiently small 
a. Then, for the function (j){x), we define (j>(9,X) by the 
relation 



IX) 



(B4) 



with [9/a\ = i and \XLjl 2 -9jl\ = i x . Here [^J is the 
Gauss notation that represents the maximum integer less 
than x. (Mathematically speaking, we should argue the 
limit a —> and a class of functions carefully, but this 
argument is not considered in this study) 

Here, neglecting the irregularity arising from the Gauss 
notation, we write conventionally 



X~ -(lii+ai ). 



(B5) 



Using Eq. I)B2|) . this implies X ~ ex. That is, the co- 
ordinate X thus defined is the large-scale coordinate de- 
scribing the long distance behavior. Next, we explain 
differentiation, integration, and noise for the functions of 
(9,X). 

differentiation: Let 4>{ x ) be a smooth function. The dif- 
ferentiation of 4>{x) is approximated by {4>i+i — </>i)/ a - 
We then have 



(B6) 
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for < io < Nq — 2. From this, we derive 

<j>(x + a) - (j){x) ~4>U + a,X+ -a) - 4>{9, X) 



d e <t>{9, X)a + T ad x 0(Q , X) + 0{a 2 ), 

(B7) 



L 



where the approximation in the first line originates from 
the discretization error and the fact that the irregularity 
of the Gauss notation was not considered. Further, in the 
second line, we have treated </>(#, X) as a differentiable 
function that might be allowed in the appropriate limit 
a — > 0. From Eq. I|B7|) . we obtain 

d x <t>{x) = {de + edxW,X). (B8) 

integration: Let (j>(x) be an integrable function. Then, 
we calculate 



dx<p(x) ~a 

N -l 

=«EE 



i=0 
N -l JVi-1 



io=0 ii—0 

Na-lNt-l , ^ 

ta ^2 <f> (ioa, —(£i! + ai ) 

i —0 ii=0 ^ 

yT-^^), (B9) 
io 1 Jo 6 

where the irregularity originating from the Gauss nota- 
tion is not considered (the third line) and the limit a — ► 
and e — > is taken (the forth line). 

coordinate dependent noise: Let £(x, t) be Gaussian white 
noise satisfying 

(£(x, t)£(x',t')} = 2S(x - x')6(t - t'). (BIO) 



We set 



a(i+l) 



dx£,(x, t) 



(Bll) 



with small a. For this discretized noise, we define (jUaM W 
and <j){6, X, t) in the same manner as that for the case in 
which <f>i — <fi(ia). Conventionally, we denote <p{9, X, t) as 
£(9,X,i) (see 1111 A"|) . Based on the definitions described 
above, we can derive 



(Z(6,x,t)t(e',x',t')) 

= 2e£S(9 - 9')8(X - X')5(t - t'). 



(B12) 



Let <p(9) be a smooth function that satisfies <p(0) = 
ip(£). We denote the Stratonovich product of ip(9) and 
£(9, X, t) as ip(9)o£(0, X, t). This product can be written 



by using the discretized form 



with an additional 



definition (pNoM = ^o.h- From this, we obtain 

d9d e [ip(9) o f (0, X, t)} = 0. (B13) 
This formula is used for obtaining Eq. (|53|l . 



APPENDIX C: PROOF OF EQS. (l68ll-(f7Tl> 

In this Appendix, we present the proofs of Eqs. 168L 
|JS3J, lf70jl. and (J7TJ- We first prove the key equality 



d ( b(9) 



e d9> 



I+(9). (CI) 



d9\ Ps (9)J yj I 
All the equations can be derived from Eq. ijClfl . 

1. proof of the key equality 

We first derive an explicit form of b(9) defined in Eq. 
(EU). Applying M (0 ^ given by Eq. (J2J to both the left- 
and right-hand sides of Eq. (|61|l . we obtain the following 
differential equation for b(9): 

d ( v s 



T dlogpM T db(9)\ 

d 9 \ Ps (ej m - -^9~ m + -^r) 



Vs(Ps(9) - 1) 



T d Ps (9) 
7 dO 



(C2) 



Integrating Eq. (|C2f> . we obtain the following first-order 
differential equation for b: 

v Tdlo gPs (0) m + TdH^ 

p s 7 d9 7 d6 

= v s [H{9)-9-K 1 ]--p a (9), (C3) 

7 

where K\ is a constant whose value is determined later. 
Here, H(9) and V(9) are defined as 



H{9)= [ d9' Ps (9'), 
Jo 

V{9) =U(9) - f9. 
We introduce b in the equation 

b(9)=p s (9) [H(9)-9 + b( 
and then rewrite Eq. (|C3|) as 



(C4) 
(C5) 

(C6) 



v s [b(9) + K x ] + -p s (9) ( Ps (9) + = 0. (C7) 



Using the equality 



v s p s eP v = -— p s ^- (p s e 0V ) , 



(C8) 



j'"d9 

we obtain the solution of Eq. (|C7|1 as 

b{9) = -K 1+ p s (9)^ V{e) [K2 - G{9)], (C9) 

where K% is a constant that is determined later, and G(9) 
is defined as 



G{9) = / dfl'e-W'). 
Jo 



(C10) 



14 



Substituting Eq. (|C9|) into Eq. IIC61). we write 



b(0) 
Ps(0) 



= H(9) -9-Kt + Ps (ey v ^[K 2 - G(9)} 



(Cll) 



Now, K\ and K2 are determined from the conditions 
(b, $ ) = and b(0) = b(£). The results are as follows: 



Ko = 



1 



1 - ePft 



G(£), 



(C12) 



d6 



(C13) 

Next, we note the identity 



Comparing Eqs. (|C16|I and (|C17|) . we obtain Eq. (|Clf> . 



2. proof of Eq. (T7CT1) 



We can rewrite D in Eq. as 



D=v, 



d9 b(0) T f e d9 



£ Ps (9) 7 Jo t 



d ( b{9) 



d9 \ Ps (9) 



Ps(0). 
(C18) 



Then, substituting Eqs. (|€TBj> and (JHU into Eq. (|€T%j) . 

we obtain 



d#V(#V /e ' - (1 - c«' f ) / d9'W)e^ 9 ' 

Jo 



e^ e / d0>(#' 
Jo 



J3fB' 



(C14) 



for an arbitrary periodic function (f>(9) with period I. Set- 
ting = e _/3c/ in Eq. I|C14|) . we obtain 



K 2 - G(0) 



1 



1 - eW 



(C15) 



Then, substituting Eq. (|C15|) into Eq. (|C11|) and multi- 
plying v s to both the left- and right-hand sides, we obtain 

Vs ^?l =v s [H{9)-9-K x ) 



T 



1 \Jo 



de 1 



J_(0') Ps (9)I + (9), (C16) 



where we have used Eq. (|15fl . 

On the other hand, we rewrite Eq. (|C3I) as 



i 



d ( b{9) 



d9 \ Ps (9) 



1 



■—b + v s (H(9) -9-Ki). 

Ps 



(C17) 



f l dO 

D=v a / — [H(9)-d]-v s K 1 . (C19) 
Jo " 

Subst ituting Eq. <|C13|) into Eq. (|C19|) and using Eq. 
'15|) . we rewrite (IC19II as 



D = -v s J o ™p 2 Me 



1 



e W - l 



-pv(e) 



(C20) 



where we have used the identity 



J ^-\p a {9)-l][9-H{9)] = 0. (C21) 



Finally, substituting Eqs. lO and (JTU) into Eq. (|C20|) . 
we obtain Eq. J7U|) . 



3. proof of Eqs. (J68J), (J5SJ|, and (17H) 



We first calculate the correlation function of rj2(X,t) 
as follows: 



( ?72 (X,^ 2 (X' ) t')) = - ^ 



T f l d9 r d9' % (9) ^ (0') 



7 Jo Uo t VpM VpMO 7 ) 

T f e d9 f e d9' $ (S)*o(fl' 



Ul 2 {9,X,t)U0\X',t')) 



i Jo t Jo * vSH 

9T 

— e 2 <S 2 (X-X')^™0, 

7 



2^e 2 (5 2 (X - X')<*(0 - 0')5(t - t') 



(C22) 



15 



where we have used (1B12|) . Using Eqs. (|64() and (|67|l . we Next, we consider the correlation function of ((X, t) 
obtain Eq. I|t)9|) . Vi(X,t) as follows: 



{[ax,t) + m(x,t)][C(x'J) + m(x'J)]) 



t r do r Ad' 



7 Jo I Jo 1 



2T r ae 
TJo ~ 



d ( 6(0) 



de \ Pb (o) 



i 



d ( b(6') 



1 



d f 6(0) 
d9 (jtfff) 



1 2 



d6> \ Ps (0' 
Ps (6)e 2 5 2 (X -X')5(t-t>). 

I 



Vps(e)p s (o') (a (e, x, i)a (^', x', o) 



(C23) 



This corresponds to Eq. 1681) . 



Because the right-hand side of Eq. IJC24I) is invariant for 




Finally, substituting Eqs. IjClfl and l|15|) into Eq. (|68|) . the exchange of i+ and J_ (see Ref. Qjj), Bn in Eq. 



-I-(0)[I + {e)f ■ (C24) 



12411 is equal to D in Eq. 1(70)1 . This corresponds to Eq. 
CD- 
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